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We develop a formalism to characterize the redshift evolution of the dark energy potential. Our 
formalism makes use of quantities similar to the Horizon-flow parameters in inflation and is gen- 
eral enough that can deal with multiscalar quintessence scenarios, exotic matter components, and 
higher order curvature corrections to General Relativity. We show how the shape of the dark energy 
potential can be recovered non parametrically using this formalism and we present approximations 
analogous to the ones relevant to slow-roll inflation. Since presently available data do not allow a 
non-parametric and exact reconstruction of the potential, we consider a general parametric descrip- 
, tion. This reconstruction can also be used in other approaches followed in the literature (e.g., the 

■ reconstruction of the redshift evolution of the dark energy equation of state w(z)). Using observa- 

f~* > ; ' tions of passively evolving galaxies and supernova data we derive constraints on the dark energy 

£\J , potential shape in the redshift range O.f < z < 1.8. Our findings show that at the la level the 

potential is consistent with being constant, although at the same level of confidence variations can- 
not be excluded with current data. We forecast constraints achievable with future data from the 
Atacama Cosmology Telescope. 



o 



PACS numbers: 



> 

Q\ • I. INTRODUCTION 

Recent observations Q, Q indicate that ~ 70% of the present-day energy density of the universe may be made 
of a dark energy component. The two leading explanations of dark energy are a cosmological constant or a slowly 
rolling scalar field e.g.,0, 0, 0, |(| but an explanation in terms of modifications to the Friedman equations(e.g.0, Q 
is also possible. In both cases this component has a negative pressure thus inducing an accelerated expansion of the 
'"^ Universe. 

A significant observational effort is directed to unveil the nature of dark energy fe.g..[50j. |5lj . H2' [53j . 
I ■ With few exceptions 0, U H El , current constraints on the nature of dark energy mostly measure the integrated value 
over time of its equation of state parameter (w = p/p) (e.g.,0, BEiEJEI) or, alternatively its energy density as a 
> ■ function of time (e.g. -JlEl Il6 | ). which depends on the integral of the equation of state parameter. These constraints 
are very tight (e.g. @j finds w = —0.98 ± 0.12) and are centered around the expected value for the cosmological 
constant, but, as pointed out by 0,0], the finding that the time average value of w is consistent with —1 does not 
exclude the possibility that w varied in time. Therefore, it is an open challenge to determine whether dark energy is 
, a cosmological constant or a rolling of a scalar field. 

A recent review of the current status of our knowledge of the observational determination of w and possible 
theoretical models to explain it is given by Q ■ 

From a theoretical point of view, it is not only important to clarify whether this energy component is dynamical or 
constant, but, in case it is not a cosmological constant, it is also of great interest to constrain the potential of the rolling 
scalar field. Since different theoretical models are typically characterised by different potentials, a reconstruction of 
the dark energy potential from observations can yield more direct constraints on physically motivated dark energy 
models. 

In this paper we present a non-parametric method to reconstruct the redshift evolution of the potential and kinetic 
ene rgy densities of the dark energy field. Our formalism introduces quantities similar to the Horizon-flow parameters 
|l9l |20| in inflation. It has the nice feature that it is easily implemented in the presence of higher order curvature 
corrections to General Relativity and different types of energy contributions in Einstein's equations, as we do in Section 
III Al Our exact reconstruction formulas determine the value of the potential at a given redshift once the matter density, 
Hubble parameter (H) and its first derivative (H) are experimentally measured at that redshift value. We discuss the 
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observational challenges to reconstruct the potential in this fully non-parametric way due to the difficulty in measuring 
H . As current data is not good enough to determine if, we present a general parameterization of the potential, based 
on an expansion in Chebyshev polynomials. In this approach, the scalar potential function at a given redshift is 
expanded in Chebyshev polynomials, which constitutes a complete orthonormal basis on a finite interval, and have 
the nice property to be the minimax approximating polynomial. Our reconstruction equation becomes a differential 
equation for the Hubble parameter, which we solve analytically, and the coefficients in the Chebyshev expansion 
become the parameters to be constrained from observations of the Hubble parameter. Our general parameterization 
can apply to other approaches that were already considered in the literature, such as expansions of the equation of 
state and we show the correspondence to some parameterizations that have been proposed in the literature. Using 
current data (in particular with recent supernovae data and relative ages of a sample of passively evolving galaxies) 
we reconstruct the potential of dark energy using our parameterization up to z ~ f .8. The reconstructed potentials 
obtained from galaxy ages and SN are consistent. Since these two data sets rely on independent physics and are 
affected by completely different systematics, this finding suggests that possible systematics are not a crucial issue. 

The reconstructed potential is consistent with being constant up to the maximum redshift of the observations, 
although current constraints do not exclude a variation as a function of redshift. We show that data obtained with 
the Atacama Cosmology Telescope will be able to greatly improve current constraints. 



II. METHOD 



A. Dynamics of the scalar field of dark energy 



The classical effective action that we shall use to describe the dynamics of the universe is 

, f ml ^ 



S= I dtd A x^{- T £{R + f(R,R^R^,...) + »—d ll qd u q-V{q)\+S somccs , (1) 



16?r v . , v , 2 

where m p stands for the four dimensional Planck mass and g^ v for the components of the four dimensional metric, 

ds 2 = dt 2 - a 2 (t)dx 2 , (2) 

which we shall consider to be an homogeneous, isotropic and spatially flat FRW cosmology, as supported by recent 
data 0. 

^sources stands for the classical action describing the physical energy content, such as matter and radiation, but it 
could also include more exotic sources (e.g. defects, cosmic strings etc.). Note also that we have implicitly assumed 
the existence of a single canonically normalised quintessence scalar field q(t, x] subject to the potential V(q). Thus, 
we have assumed that this potential is independent of the derivatives of the scalar field. For generality, in Eq. (JTJ 
we include the effect of higher derivative terms in the gravitational sector of the theory p| • These are described by 
the function f(R, R^R^, . . . ) of the different invariants that we can construct out of the metric and its derivatives. 
In four dimensions, the most general lowest order corrections to Einstein's classical action would be described by 
/ = f3 R 2 + SRVR^v |5(| (see, for example, |2*Hl. Other corrections that have been considered, include arbitrary 
functions of the scalar curvature f(R) |22], which include as particular examples linear combinations of negative 
powers of these invariants [23| . 

We focus on cosmologies given by (J2J, and shall restrict ourselves to classical configurations q = q(t), configurations 
that do not break the homogeneity and isotropy of spacetime. The energy momentum tensor of this scalar field 
configuration is that of a perfect fluid, with density p q and pressure p q given by 

Pq = K(q) + V(q) , Pq = K{q) - V(q) and K = l -q 2 . (3) 

where K denotes the kinetic energy of the field. Under these assumptions, one is led to consider Einstein's equations, 
plus the Klein-Gordon equation of motion for the scalar field. The first ones reduce to Friedmann's equations 

H 2 = | (PT + P q ) , 

- = - g (PT + iPT + Pq+ 3Pq) , 

where k — 8Tr/m 2 (or n = 8irG). In Eq. J3J we introduced the compact notation and pt for the total energy density 
and pressure. For example px denotes the full energy density contribution of Sources and of the higher derivative 
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curvature terms f(R, R^ V R^ V , . . .). Thus if the sources are a collection of n perfect fluids with constant equation of 
state u>i i = 1, . . . ,n, px and pr are 

n n 

PT = ^Pi + Pf, Pt ^^2^iPi+Pf , (5) 

i=l i=l 

where and pf describe the contribution from the higher derivative curvature terms. For the particular function 
f(R 2 , R^Rpv) introduced above, such terms would be written as 

(3) / •* \ 2 

Kpt = -12H — (3(3 + 6) -6 ( - ) (3/3 - 5) + 6ff 4 (15/3 + 75) - 36ff 2 - (J3 + 5) 
a \a J a 

a /a\ 2 

np f = 4 ((5 + 3/3) + 12# (5 + 2/3) + 8ff 2 - (6 + 9/3) + 2 - (6+ 3/3) + 2# 4 (15/3 + 75) . 

a a a \ay 

On the other hand, the scalar field q(t) equation of motion reduces to 

q + 3Hq + V' = 0, (6) 

where V' — dV/dq. 

B. Reconstruction procedure 

In this section we provide exact analytical expressions in which both, the kinetic and potential energies of the 
quintessence field q(t) 7 depend on quantities more directly observable such as the energy densities, the Hubble constant 
H and its derivatives. Although the higher curvature corrections are not directly observable, they will also have to 
appear in the expressions: they can be taken into account for a given model that is for a given parameterization 
of the functional /. Provided one has an independent way of determining the densities, H and H, the value of the 
potential V(z) at a given redshift z where these measurements are available, can then be fixed, up to experimental 
uncertainties. If also higher order derivatives of H are known, higher derivatives of the potential can be determined 
rig oo wn i cn can he used to probe the flatness of the potential. 

We use the analogous of the inflationary horizon-flow parameters 0,|2(j {£ n }, which are defined recursively by 

dlog | e n | 
e« + i = —w— , n>0 

where N — log(a(t)/a(t i )) is the number of e-foldings since some initial time i, and Eq = H(Ni)/H(N). There are 
many similarities between the period of inflation and the present-day accelerated expansion, but, despite the fact that 
inflation happened 13.7 billion years ago, and the accelerated expansion is happening today, as we will see, it is not 
observationally easier to reconstruct the dark energy potential than it is to reconstruct the inflationary potential. In 
the equation describing inflationary dynamics the contribution due to matter can be ignored, but it can't be ignored 
when describing today's expansion. Moreover, the detailed shape of the primordial power spectum from CMB scales 
to large scale structure scales, and the nature of the primordial perturbations offer a window to test the last 4 inflation 
efoldings; conversely, in the case of dark energy, dark energy started dominating at z < 1, and between then and 
now the Universe expanded only by a factor < 2. In addition we can measure with exquisite precision perturbations 
from inflation but have not detected perturbations from dark energy, which is a very challenging task [24j]. On the 
other hand we do not have strong constraints on the energy scale of inflation, that is on the "normalization" of 
the inflationary potential [25l |2(j but as we will see, since the matter content of the Universe can be independently 
determined, for a flat Universe, we have some constraints on the quintessence potential normalization. 

Keeping in mind the different kind of challenges that a quintessence potential reconstruction faces, we proceed with 
our program. For our purposes, it will be useful to have explicit expressions for the first two parameters 

H d _dH (l + z) 
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which, we will show, are needed to determine V and V 57]. We use the second Friedmann equation {3} to express 
the first horizon-flow parameter E\ in terms of the energy and pressure densities: 

3 p T + Pq + PT + Pq 
2 PT + Pq 

If we write {p q , p q } in terms of its kinetic and potential energy components, as in @, we can use (JHJ to express, 
e.g., the kinetic energy in terms of the potential energy as 

i.2 i r , . •■-! 



2<? = — — |ei (pt + VO - 2 (Pt + pr) j . (10) 

Finally we can use the first Friedmann equation 0) to solve for the value of the kinetic energy and the potential at 
a given redshift z 

K(z) = \e=e x ^- 1 -{p T +p T ) , (11) 



V(z) = (3- £l ) — + Up T -p T ). (12) 

K A 

This is the generalisation of eq. (16) in [§o| which was derived in the context of inflation. 

Equation (|12|l is a general and exact reconstruction formula for the potential of a quintessence field given the 
assumptions followed in this paper. 

Here, we shall focus in the constraints on the potential at redshifts smaller than 1000. Therefore, we shall neglect 
the radiation energy density contribution. Furthermore, if we also neglect the contribution from the higher-order 
curvature terms, the expression for the potential simplifies 

V{z) = {i-e 1 ) — -\p m . (13) 
Analogously, for the kinetic energy we obtain 

K(z) =ex — -\ Pm . (14) 

K Z 

Ideally, our goal would be to constrain the functional form of the potential, V[q], and this is not what (|12fl provides. 
If the function q(z) was known this would be straightforward, but q{z) is not an observable quantity. We will show 
later that V[q] can be obtained if equation (|14fl can be integrated. 

We can next determine the first derivative of the quintessence potential. We rewrite 10 as 

V =-{q)- l {ZHq 2 + qq} . (15) 

where all terms are already known, except for q q which can be obtained from the time derivative of the kinetic energy 
(|ll|l . The end result can be expressed as 



4^ 1 2H 2 £l ^ y 1 1 6 3 6# 3 ei 



3H (p T + pt) + x [Pr + Pt) 



(16) 



Thus, if the values of Pt{z),pt{z) and {H, £i, 62} or equivalently, {H, H, H}, can be experimentally determined for 
some redshift z, (|16JI yields the first derivative of the potential V'(z). As in the previous discussion, the determination 
of V would require the knowledge of q(z). 

The above formula is the exact result given some energy density content pt, with associated pressure px- If wc 
restrict ourselves to a single matter component and neglect the higher order curvature terms, the first potential 
derivative reduces to 

V " ,3) = - 3 ^ ,,2(E ' ) " 2 { 1 -5lfe^}" /2 { 1 + ?-T-Ilfe'""} ■ (17) 

and l(T5jl reproduces equation (17) in [2(| when the matter density vanishes (p m = 0). 

In this case, the first derivative of the potential V'(z) is known if one can measure p m (z = 0), {H, H, H}. 

Analogously exact expressions for higher order derivatives of the potential S r ^V[q]/dq r can be obtained by taking 
the time derivative of (|16fl and using the exact expression for the kinetic energy (|11|) . 
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C. Redshift parameterisation of the potential 

In section Til Bl we have shown that an exact reconstrution of V(z) is possible only if H(z) and H(z) are known. 
While the determination H(z) is an observationally challenging task (e.g, [13, and llVJI . the determination oiH(z) 
is even more formidable. In this section we shall not attempt a non-parametric and exact reconstruction of V(z), we 
shall instead consider a parametric description of the potential (V(ai, z)) in terms of the redshift z and parameters 
cti. In section IW1 we will then use currently available observations to constrain the potential parameters and discuss 
future prospects. Hereafter we will set pf = and defer the more general case of pf ^ to future work. 

Equation {T3J, can be rewritten in terms of the independent variable z as 

3H\z) -\{l + z) = k (v(oi, z) + lp m (z)) = g(oi, z) . (18) 



(19) 



2 v ' dz \ v ' ' 2 

This is a first order non-linear differential equation which can be integrated analytically: 

H 2 {cti, z) = Hl(\ + zf - 2{l + zf [ gfa, x) (1 + x)- 7 dx 

Jo 

= (ff 2 -^p m ,o) (1 + z) 6 + ^p m (z)- 2(1 + z) 6 V( ai ,x)(l + x)- 7 dx . 

Hereafter the subscript denotes the quantity evaluated at z = 0. 

In this approach if we now consider the kinetic energy of the quintessence field we obtain a first-order non-linear 
differential equation for q(z) 

i • r z 

{qf = (1 + z) e V - 6(1 + zf / V{a h z)(l + z)- J dz + K (l + zf . (20) 



2 



o 



or equivalently 



1 (dq s 2 



(1 + z) 1 H 2 {a u z) = i^ 1 H 2 (a t , z) - p m (z) ~ V(a h z) , (21) 



2 \dz 

which can be integrated to obtain q(z) and thus V^[o!i, q] from V(oti, z) 



q(z) - q(0) = ± [ Z — dZ {6k- 1 H 2 (a u z) - 2p m (z) - 2V(a t , z)} 1/2 , (22) 

J (1 + z) H(ai, z) 

where the ambiguity in sign comes from the quadratic expression for the kinetic energy. Typically if we think of an 
scalar field rolling slowly along its potential, the plus sign will be the relevant one. 
For example let's consider a simple two-parameters parameterization of the potential: 

V = \(l + z) a , (23) 

which yields 

H 2 {z) =H (l + z) 6 ~2I a (l + z) 6 , (24) 

where 

la = ™Pm(Q) [(1 + z)- 3 - 1] + [(1 + z)"- 6 - 1] , a ? 6 (25) 

o a — b 

h = ~Pm(Q) [(l + ^r 3 -l] +A/elog(l + z) , a = 6 (26) 

If we can neglect the kinetic energy (that is if a << 1) then this potential correspond to a constant equation of 
state w: a = 3(w + 1) as in the Ratra- Peebles Q case and f2 9j o = 6A/[p c (6 — a)]. 
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D. Chebyshev reconstruction 

An interesting parameterization of the potential involves the Chebyshev polynomials, which form a complete set of 
orthonormal functions on the interval [—1, 1] . They also have the interesting property to be the minimax approximating 
polynomial, that is, the approximating polynomial which has the smallest maximum deviation from the true function 
at any given order. We can thus approximate a generic V(z) as 

N 



V(z) ~ Y, X nT n (x) (27) 



71=0 



where T n denotes the Chebyshev polynomial of order n and we have normalized the redshift interval so that x — 
2z/z max — 1; z max is the maximum redshift at which observations are available and thus x G [—1, 1]. Since |T„(a;)| < 1 
for all n, for most applications, an estimate of the error introduced by this approximation is given by Ajv+i- With 
this parameterization, the relevant integral in i|2U|) becomes: 

V{y){l + y)- 7 dy=^-Y\ n T n (x){a + bx)-' ' dx = Va„F„(z) (28) 

J ° Z n=0 J-l n=l 

where a = 1 + z max /2 and b = z max /2. These integrals can be solved analytically for any order n as shown in the 
appendix: F n are known analytic functions which are reported in appendix A. 
We obtain: 



N K ( i 



(29) 



where p c denotes the present-day critical density. 

Equation l|29l) seems to describe the potential with N + 1 parameters (Ao ■ ■ ■ Ajy). However, since we assume a flat 
Universe, fl m ,o constrains the coefficients of the Chebyshev polynomials in expansion 127|) and the kinetic energy of 
the field. For example, if the potential is constant (i.e in a cosmological constant case) then it is completely described 
by only one parameter, Ao, and since q = 0, Vo + Kq = Vq = p q ,o = Q q ,opc, we have that Ao = f2 g ,oPc- However if the 
potential is not constant we obtain 

\q 2 = V Q (l + zf-V(z)-6{l + z) 6 f V(z)(l + z)- 7 dz + K (l + z) 6 (30) 

6 







2 

= V (l + zf - V{z) - 3(1 + zfz max ^nF n + K (l + z) 



and since Vq — Y^,iLo At(— 1)™ we have the constraint: 



5^M-i) n + ^o = n,,oPc. (3i) 

i=0 

In section [Tvl we show the constraints that can be obtained for the first few Chebyshev coefficients from currently 
available data. 

III. EQUATION OF STATE 

It is widespread to parameterize dark energy not by the scalar field potential but by its equation of state. In this 
section we connect the two descriptions. 

A. Time dependent density scaling 



Standard contributions to the energy momentum tensor in Einstein's equations are characterized by a parameter 
that governs how their energy densities decrease with the expansion of the universe. For the energy density of the 
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scalar field we can write in all generality 



«o 



i{t) 



(32) 



The Klein-Gordon equation 10 can be expressed as the conservation equation for the energy momentum tensor 
describing the scalar field, 

p q + 3H (p q +p q ) =0 . (33) 
Thus, the pressure of the scalar field p q (t) can be expressed as a function of the time-dependent exponent ^(t) 



lit) 



1 1 ( a ° \ TJ-l 



p q (t) 



Using ansatz l)32p. we obtain the kinetic and potential energies for the scalar field 



(34) 



p q (t) At 



Q ) 



V 



Pq [1--Aw q (t) 



where we introduced the function (compare to |29^ 



Aw q (t) 



7(t) - H- 1 log 





c?7 


(a(t)) 





which depens on the ratio between kinetic and potential energies of the scalar field via 

2 



Aw q (t) = 



1 + V[q]/K[q] ■ 



(35) 
(36) 

(37) 
(38) 



The function Aw q (t) controls the deviations from the equation of state of the scalar field w q (t) from being exactly 
-1 i.e. a cosmological constant. Indeed, from p q (t) = w q (t) p q (t), we have 



w q {t) = -1 + Aw q (t) . 



(39) 



Thus constraining Aw q (t) is equivalent to constraining the time (or redshift) evolution of the dark energy equation 
of state for which there are different independent observational constraints (e.g., [3ll30|. and references therein) 
We can relate ^(t) to Aw q (t) by integrating l|37|) . Without loosing generality, we do this in terms of the redshift 



70) 



log( 



dy 



-, r / Aw n (v) 

(l + z)J qyy, l + y 



(40) 



where we did not include a possible integration constant, since it is physically irrelevant, i.e. it just redefines the value 
of the energy density of the scalar field today (p q (Q)). This relation provides an expression for the potential energy 
density of the scalar field once the equation of state is known, using iffifl) : 



V(z) = p g (0)e 3 ^ Aw * iy) t& (l - ^Aw q (z)j 



(41) 



The case of a constant Aw q , or equivalently, a constant equation of state w q would correspond to the redshift 
parameterization considered in 1)23(1 with 



and a — 3(1 + w q ) 



X - Pq (0) [ 1 - -Aw q 



Thus, for a physical situation resembling a cosmological constant (w q ~ —1), the parameter a <JC 1, as we claimed in 
the previous section. Some constant scalar field energy density scalings were discussed by |4|. The above formalism 
is a natural generalisation of their models. 
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B. Beyond a rolling scalar field 



Up to now we have assumed that dark energy is given by a scalar field. The techniques developed here however 
can also be used in a more general context. 

If the dark energy is not due to a scalar field we can still describe it as a fluid with a given equation of state. As 
long as the equation of state w is > — 1 the two descriptions are equivalent. For example for a given set of parameters 
Aj and H and £l m .o, V(Af, z) is known, the kinetic energy can be computed using eq. 1311 thus Aw is known (eq. I38|) 
and so w(z) (eq. EHJ). 

However if we want to allow w < — 1 then the scalar field description as presented here fails. However, analogously to 
section lil Dl we can expand the redshift dependence of w in Chebyshev polynomials without imposing any restrictions 
on the values that w can take. We thus obtain 



N 



i=0 



and 



H 2 (lu. 1 ,z)~H%(1 + z) 3 



N 



Sl mfi + Q qfi exp ^2 u n G n (z) 



n=0 



(42) 



(43) 



where G n is the analogous of F n of section III Dl F n is a linear combination of integrals Jj, i = {0, n} while G n is 
the same linear combination of the integrals Ji , i = {0,n} (see appendix). 
Note that in this parameterization the present-day value of w is given by 



N 



= £(-d 



(44) 



i=0 



Given the parameterization (|42|l of w(z) subject to the constraint w > — 1 one can always obtain: 



V{z) = p qfi {l + z )W+<*>) 



exp 



N 



5 ' tOiGi(z 



/V 



1 -^2<jJiTi(x(z)) 



i=0 



(45) 



In section Hvl we show how currently available data can be used to constrain the first few Chebyshev coefficients of 
this expansion. 

In the remaining of this section we will compare some models presented in the literature with the parameterization 
presented here. Clearly, the case of a constant equation of state corresponds to u>.- L — for % > 0. The linear 
parameterization in z pH l32| corresponds to u>i = for i > 1, and in particular wo = ^0 — uj\ and w' — 2uj±/ z max . 
Finally the linear parameterization in a I33l 1341 . w — wq + w a z/ (1 + z) for w a <C wo can be closely approximated by 
LOi = for i > 2, with the constraint 144f) . [35| pointed out that a simple, 2-parameter fit may introduce biases: the 
expansion (|42ll allows one to include more parameters by increasing N as the observational data improve. 



IV. OBSERVATIONAL DETERMINATION OF H(z) 

Section lrmi has illustrated that it is necessary to determine observationally H(z), H(z),H(z) in order to reconstruct 
V[q] and its first derivative V'[q\. Here we present a determination of H(z) based on the method developed by j3|J 
and we emphasize the difficulties of computing H(z), H(z). We also present the constraints that can be achieved on 
the evolution of the quintessence potential and the dark energy equation of state from present and future data. 



A. Differential ages of passively evolving galaxies 



The Hubble parameter depends on the differential age of the universe as a function of redshift in the form 



(46) 



9 




250 



200 



150 



1.0 1.5 2.0 




100 



1.0 1.5 2.0 



FIG. 1: Left panel: the absolute age for the 32 passively evolving galaxies in our catalogue (see text for more detals) determined 
from fitting stellar population models is plotted as a function of redshift. Note that there is a clear age-redshift relation: the 
lower the redshift the older the galaxies. Right panel: the value of the Hubble parameter as a function of redshift as derived 
from the differential ages of galaxies in the left panel. The determination at z ~ 0.1 indicated by the '+' symbol is the hubble 
constant determination of H from |2?j |. The dotted line is the value of H(z) for the LCDM model. 



Therefore a determination of dz/dt directly measures H(z). In [22| we demonstrated the feasibility of the method 
by applying it to a z ~ sample. In particular, we used the Sloan Digital Sky Survey to determine H(0) and showed 
that its value is in good agreement with other independent methods (see |27J for more details). With the availability 
of new galaxy surveys it becomes possible to determine H{z) at z > 0. Here we use the new publicly released GDDS 
survey ^3 an d archival data [H, SES IMIi2. E3 to determine H(z) in the redshift range 0.1 < z < 1.8. We proceed 
as follows: first we select galaxy samples of passively evolving galaxies with high-quality spectroscopy. Second, we 
use synthetic stellar population models to constrain the age of the oldest stars in the galaxy (after marginalising over 
the metallicity and star formation history), in similar fashion as is done in |2?| |. We compute differential ages and use 
them as our estimator for dz/dt, which in turn gives H(z). 

The first sample is composed of field early-type galaxies from j^l, |3^, |4jj . In |53 we derived ages for this sample 
using the SPEED models 0. The second sample is from the publicly released Gemini Deep Survey (GDDS) 37]. 
GDDS has high-quality spectroscopy of red galaxies, some of which show stellar absorption features, indicating an old 
stellar population. The GDDS collaboration has determined ages (and the star formation history) for these galaxies 
|45| : they conclude that for a sub-sample of 20 red galaxies the most likely star formation history is that of a single 
burst of star formation of duration less than 0.1 Gyr (in most cases the duration of the burst is consistent with Gyr, 
i.e. the galaxies have been evolving passively since their initial burst of star formation). To determine the galaxies 
ages they use a set of stellar population models different than SPEED. We have re-analize the GDDS old sample 
using SPEED models and obtained ages within 0.1 Gyr of the GDDS collaboration estimate. This indicates that 
systematics are not a serious source of error for these high-redshift galaxies. We complete our data set by adding the 
two radio galaxies 53W091 and 53W069 |4l|, |42, Hj] . In total we have 32 galaxies. 

Fig. n (left panel) shows the estimated absolute ages for galaxies in the above samples and their ler error bars. 
There is a distinguishable "red envelope" : galaxies are older at lower redshifts. 

The next step is to compute differential ages at different redshifts from this sample. To do so we proceed as follows: 
first we group together all galaxies that are within Az — 0.03 of each other. This gives an estimate of the age of the 
universe at a given redshift with as many galaxies as possible. The interval in redshift is small to avoid incorporating 
galaxies that have already evolved in age, but large enough for our sparse sample to have more than one galaxy 
in most of the bins. We then compute age differences only for those bins in redshift that are separated more than 
Az = 0.1 but no more than Az — 0.15. The first limit is imposed so that the age evolution between the two bins is 
larger than the error in the age determination. This provides with a robust determination of dz/dt. We note here that 
differential ages are less sensitive to systematics errors than absolute ages (see |44| for detailed discussion, specially 
their table 2). The value of H(z) is then directly computed by using Eq. 0^1 This is shown in fig. ^ with ltr error 
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FIG. 2: Regions in the \\/p c vs Ao/p c (left panel) and A2/ p c vs Xo/pc (right panel) excluded at the 1 — a and 2 — a joint 
confidence level, by the priors and the constraints that the kinetic energy in the quintessence field must be positive and that 
at all redshifts p m + p q must be positive. 



bars. Also shown (dotted line) is H{z) for the LCDM model. 



B. Constraints on the potential 

Following the discussion in sec III Dl we present constraints on the shape of the potential achieavable from present 
and future data sets. Figure J3J shows the constraints on the first three Chebyshev coefficient for the potential that can 
be obtained from our galaxy sample, combined with the determination of the Hubble constant at z = 0.09 obtained 
by from the SDSS luminous red galaxies. We have assumed a flat Universe and marginalized over a gaussian 
prior on S} mfi (Q m ,o = 0.27 ± 0.07 (e.g., 46]) and a flat prior on H (30 < H < 100 Km/s/Mpc). We have used 
only the large scale structure prior on f2 mj o, as the determination of [46j is insensitive to dark energy. Conversely, 
CMB constraints on the matter density of the Universe are highly sensitive to the assumptions about the nature of 
dark energy (see e.g., in particular figure 12), and thus should not be used in this context. Of course, the addition 
of CMB data can greatly improve the constraints on the nature of dark energy, but this need to be done in a joint 
analysis and it is left to future work. 

Some regions of the parameter space are unphysical as they would yield a negative kinetic energy or p m + p q < 0; 
the combined effect of these priors in the Xi/p c vs Xo/p c plane and \ij p c vs Xo/p c plane is shown in fig. We 
consider only the region < \o/p c < 1.1, —0.5 < Xi/p c < 0.5 and —0.5 < A2//OC < 0.5. 

In fig. |3|we show the one and 2 sigma joint confidence contours in the \q/Pc vs X\/p c and Xo/p c vs \ij p c planes, 
obtained from our H(z) determination. When adding the HST key project prior on Hq |47| the contours remain 
virtually unchanged. For comparison in figure|Ijwe show the constraint obtained by using the recent supernovae data 

of H- 

Figure shows our best fit reconstructed V(z) from our H(z) determination (left panel) and from the SN data 
(right panel), and the 68% and 95% confidence regions. The present constraints are consistent at the 1-a level with 
a constant potential (that is the cosmological constant scenario) . 

The two determinations (one based on relative galaxy ages and one SN data) are consistent with each other. The 
two methods are completely independent and are based on different underlying physics, different assumptions and 
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FIG. 3: Constraints in the Ai/p c vs \o/p c (left panel) and \il pc vs Xo/p c (righ panel) obtained from H(z) measurement based 
on relative galaxy ages. Contour levels are 1— a marginalized, 1— a joint and 2— a joint. The diamon shows the location of the 
maximum of the marginalized likelihood. 
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A o/Pcrit A o/Pcrit 



FIG. 4: One and two sigma joint constraints in the Xi/p c vs Ao/p c plane and \2/pc vs \o/pc obtained from the Riess et al. 
(2004) supernovae data. 
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0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 



z z 

FIG. 5: Reconstituted V(z) from relative galaxy ages (left) and from Supernovae (right). The gray regions represent the 1- and 
2- a confidence regions. In the left panel the dotted line shows the constraint imposed by the prior. 

affected by systematics of completely different nature. The fact that they agree indicates that possible systematics 
are smaller than the statistical errors. 

With current data there is a degeneracy between the first two coefficients, but we can place an upper limit to the 
kinetic energy in the quintessence field today: the contribution of the kinetic term to p q is less than 40% at the 2-tr 
level and the best fit value is at 0. 

The Atacama Cosmology Telescope (ACT; 48] www.hep.upenn.edu/act) will identify, through their Sunyaev- 
Zeldovich signature in the cosmic microwave background, all galaxy clusters with masses > 10 Mq in a patch 
of the sky of angular size 100 square degrees. Thus ACT will yield > 500 galaxy clusters in the redshift range 
0.1 < z < 1.5 . For all these clusters, spectra of the brightest galaxies in the cluster will be obtained by South African 
and Chilean telescopes. This will provide us with an unbiased sample of > 2000 passively evolving galaxies from 
z = 1.5 to the present day. To estimate the performance of ACT galaxies at reconstructing the dark energy potential, 
we have estimated that we will have 2000 galaxies for which age s have been determined with ~ 10% accuracy and 
therefore ~ 1000 determinations of h with ~ 15% error. Our|58j forecasts in the reconstruction of the dark energy 
potential are shown in fig. El We have marginalized over a flat prior on the Hubble constant 30 < Hq < 90 km s _1 
Mpc -1 and a gaussian prior on fl m ,o, ft m ,o = 0.27 ± 0.035, as an estimate of the improvement of this determination 
from galaxy surveys. 

C. Constraints on the equation of state 

It is illustrative to work out the consequences of the constraints found on Ao, Ai,A2- Let's consider the potential 
(|23|l that gives rise to a constant equation of state. If a is small then one can approximate the potential with A(l + az) 
and thus identify the coefficients in the Chebyshev expansion: Ao — ► A and Ai — ► Xa, A2 = 0. We thus obtain 
A1/A0 < 0.3 at the 1-a level; since A1/A0 — a = 3(w + 1) we obtain w < —0.9 at the 1— a level. 

As illustrated in sec. IIII Bl for more general cases we can expand the redshift evolution of the equation of state 
parameter in terms of Chebyshev polynomials. Here we show how constraints on w(z) obtained from our galaxy- 
sample with the differential ages method compare with other constraints. For example in figure[7](left panel) we show 
the constraints in the plane u>o vs u>\ (i.e. we impose N = 1 in I42fl . where we have used the HST key prior for Hq 
and the prior f2 mj o = 0.27 ± 0.04 as in m. The contours show the la marginalized, la and 2a joint confidence levels. 
To compare with the SN constraints of [l| recall that their wo is ojq — uj\ . Thus the degeneracy seen in the figure is a 
constraint on wq. Three points in the luq y s toi parameter space are indicated by the diamond, star and '+' sign. These 
points are at the la joint confidence level, well within the la marginalized level and at the la marginalized level, 
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FIG. 6: Predicted constraints for a experiment with 2000 galaxies for which ages are measured with an accuracy of ~ 10%. 
The constraint in the Chebyshev coefficients (left panels; circles show the location of the maximum marginalized lielihood while 
'+' show the location of the maximum of the joint -5D- likelihood) and the reconstructed dark energy potential (right panel) 
are significantly better than current constraints (see text). We have a LCDM model as fiducial. The Atacama Cosmology 
Telescope will identify about 500 galaxy clusters in the redshift range 0.1 < z < 1.5, for at least 2000 galaxies there will be 
spectroscopic follow up and therefore galaxy ages can be derived. 




FIG. 7: Left panel: Constraints in the cjo vs ui obtained from our galaxy sample with the differential ages method. The 
contours show the 1 — a marginalized, 1 — a and 2 — a joint confidence levels. The degeneracy constrains wo = coo — wi. Three 
points in the parameter space are selected. Right panel: difference between the hubble constant in a given model and the 
Hubble constant in the LCDM model. The points with error-bars are our data points, the long-dashed line corresponds to the 
LCDM model (*-point in the left panel), the dot-dashed line corresponds to the "diamond" -point and the dot-dot- dot-dashed 
one to the '+'-point. 



respectively. In particular the '*' point correspond to the LCDM model. In the right panel we show the difference 
between the Hubble parameter for a given model and that in the LCDM case. Also our determinations of H(z) 
are shown. The long-dashed line corresponds to the LCDM case (* point), the dot-dashed line corresponds to the 
"diamond" -point and the dot-dot-dot-dashcd line to the '+'-point. It is clear that more data-points in the redshift 
range around z = 0.7 would help in breaking the degeneracy. 



V. SLOW-ROLL DARK ENERGY 



The constraints derived from our observational determination of H(z) combined with our theoretical analysis suggest 
that observations in the redshift range 0.1 < z < 1.8 are consistent, at the la level, with a cosmological constant 
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equation of state (w = —1). 

This suggests to analyse more closely the conditions under which a quintessence field could resemble such an 
equation of state in that redshift range, because this is the challenge we will be facing in the near future. 

There are at least two different approaches that one can attempt: either work with a generic potential and determine 
the properties it has to satisfy to resemble a cosmological constant, or attempt to argue some universality in the 
functional form of the potential due to its expected flatness in field space. 



A. Slow-roll in redshift 



Given a generic potential scalar field in the presence of a non-negligible matter energy density p m {z), we would 
expect the conditions the potential has to satisfy to be a natural generalisation of the slow-roll conditions during 
inflation, including the effects of matter. These two conditions are : 

«,,(*) ^^«0 Vze[0, z }. (47) 

The first one ensures that dark energy behaves approximately as a cosmolgical constant at a given redshift z, whereas 
the second ensures that such property is maintained in time. 

There are several equivalent ways of studying the consequences of these conditions. In terms of the kinetic scalar 
field energy K[q] and its potential energy V[q], Eq. I|47|) implies that V[q] ^ K [q] and that the ratio, K [q]/V[q] is 
nearly constant in time: 

«,(*)«-! => ™«1, 

V m (48) 
MM w Q ^ dK[q]/V[q] ^ 
dt ~ dt ~ 

In terms of the fundamental degrees of freedom, q(t), conditions 147(1 are equivalent to 

u> B (z)«-l l -q 2 ^V[q], (49) 

* m"^* 1 ' (50) 

where the last inequality is derived from the identity 

Under these circumstances, the first Friedmann equation and the Klein-Gordon equation © reduce to 

3^L^ Pm + v, (52) 

3Hq w -V , 

which are the extension of the slow-roll equations used in inflation in the presence of matter. One can now rewrite 
conditions l|49l) and (|50|l respectively as 

m p V'\ ,„ . 0„ 



v- y <<48 »( 1 + "!f I 

™^«24»(l + |^), (54, 

where we already used the fact that PmV^ 1 ~ flmSJ^" 1 whenever l|49() is satisfied. 

Following the discussion in section lil (^1 it is also convenient to rewrite these conditions in terms of redshift derivatives 
of the potential V[g(z)]. The analogue of conditions i|53|> and are: 

1 dV 6 tKK\ 

y l + z 

fdvy 1 d 2 V 5 
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B. Slow-roll in the field 



Phcnomcnologically, there are many inequivalent functionals that could be chosen to describe the quintessence field 
dynamics. Each of them, would typically depend on a set of undetermined parameters, which would be determined 
by fitting them to observations, as we did in Section llVl In order for a generic potential to look indistingishable from 
a cosmological constant, these parameters need to be highly fine tuned. 

It is precisely this fine tuning that suggests that, independently of the functional form of the potential, the potential 
will allow an expansion in terms of the variation of the unobservable scalar field variation A q (t) — q(t) ~ q(0), 
measuring its variation from its current value today. 

Let's assume that there is a certain period of physical time around today, i.e. t — 0, and consistent with the range 
of redshift covered in this work, where the variations in the scalar potential are small in field space. In other words, 
the potential is "flat". Under these conditions, and independently of its functional form, the potential V[q] can be 
approximated by 

V[q] « V[q(0)] + V'[q(0)}A q (t) + \ V"W)\ (M*)) 2 + Q ((A ? W)) 3 ) . (57) 

Let us emphasize that such an expansion is always viable for small enough A q , but the Taylor expansion can have a 
wider validity if the potential is flat enough, that is if the derivatives of the potential are small \ V n \ << Vb and if the 
kinetic energy is small K << Vq. 

For this to be a good approximation the following two conditions should be satisfied: 

~ K (q(t) - q(0j) « V ' , (58) 
V{ (q(t) - q(0)) < V , (59) 

where we introduced the notation V Q {n) = £^Si.[q(0)]. 

We shall also assume that the energy of the scalar field q(t) is dominated by the potential energy 



\q 2 « v[ q ] . (60) 

so that the scalar field dynamics can resemble a cosmological constant (see (1481) ) and the rolling due to the kinetic 
energy is small. In the following, we shall proceed to attempt to integrate the system perturbatively. At zeroth-order 
check in the potential expansion the first of Friedmann's equation in reduces to 

H 2 ~- [pt,o — +V 



3 \ ' \ a 

where a(0) = 1. If pr,o = p m ,o, w ~ and this reduces to the case of a LCDM universe. 
The exact solution 

iH = ^- y° + t ^\ ct , (6i) 
1 + w 1 + yo tanh ct 

and 

{a(t)f 1+w) = (cosh ct) 2 (1 + y tanh ct) 2 , (62) 
becomes, for matter + dark energy universe 

a(t) 3 = {Vo - 1) sinh2 ( ct + 5 ) , (63) 

which yields a Hubble parameter, at zeroth order H^)^ 1 = 3 tanh(c£ + c)/(2c), where we have introduced two 
dimensionless parameters: 



Ipo 

yo = y 

and the dimensional one 



yo = \/yr + l, tanhc = y 1 , (64) 



V^Vo^^ ■ ( 65 ) 
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In a matter + dark energy universe (w = 0) and in our current approximation c ~ 3HoWi\fi/2 ~ 0.09 Gyr 1 , 

whereas yo ~ ^/^m.o/^g.o + 1 ~ 1-18 and c ~ 1.2. Thus the age of the Universe is ~ 13 Gyrs. 

We can then proceed to integrate the Klein-Gordon equation by taking the first non-trivial contribution coming 
from the expansion l|57|) . and plugging in the zeroth order Hubble parameter 

q + 3Hq + V£~Q . (66) 

The solution for the velocity of the scalar filed is 

. = k - F(t) 

q (coshrf + y sinhci) 2 /( 1+,i ') 1 ' 

where 

F(t) = Vq [coshcx + y sinhca;] 2/(1+ '" ) dx = ( if w = 0) ^ F ' (y$ - 1) Q sinh 2(ct + c) - (ct + c)\ . (68) 
Here, k is the integration constant and we have used the fact that y$ — 1 > and the identity 



cosh ct + y sinh ct = y y$ — 1 sinh {ct + c) . (69) 

Note that we could identify the constant of integration k with the expression k — q(0) + F(0) involving the kinetic 
energy of the scalar field today. A second further integration yields the dynamical evolution for the field, 

q(t) - t?(0) = fci + k [ (a(x))- 3 dx- [ F(x) {a{x))~ 3 dx . (70) 
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, k / tanh c \ V Q ' ( c ct + c ■ 



which in a matter + dark energy universe becomes 

U- , 

c (uq — 1) tanh c \ tanh (ct + c) / 2c 2 \tanh c tanh (c< + c 
where 

fc = i(0)+F(0) , F(0) = {yl - 1) ^ (™^c _ ^ _ (?2) 

The approximations presented here are therefore valid if VJj'/Vb — <?(0)] < 1. Figure [S] shows, for some choices 
of V0, Vq and Kq, the range in lookback time where this approximation is valid. 



VI. CONCLUSIONS 



We have proposed to constrain the nature of dark energy, a rolling scalar field or a cosmological constant, by 
reconstructing its potential as a function of redshift. We have presented a formalism, similar to the horizon- flow 
parameters in inflation, to relate quantities characterizing the dark energy dynamics, i.e. potential and kinetic energy 
densities, to direct observables such as the matter density p m (z), the Hubble parameter H(z) and its derivatives. 
This is the core of our reconstruction programme, and our results are summarised by Q12(l. 1)11(1 and 1(16(1 . which 
provide the value of the potential and kinetic energy densities, and the first derivative of the potential, as a function 
of redshift. These expressions are valid even in the presence of higher order curvature corrections to General Relativity 
and exotic matter sources. In principle, integrating the exact reconstruction formula for the kinetic energy, allows 
one to determine the function q{z). Using the latter, one can infer the real shape of the potential V[q] from the 
determination of V(z). 

We have then focused on the case of a expanding Universe with only matter and dark energy components and at 
z <C 1000. In this case the above expressions simplify to ((13(1 . 1(14(1 and ((17(1 . This exact reconstruction formulas 
are currently difficult to be used due to the experimental challenges in determining the derivatives of the Hubble 
parameter IIVI However, given a parameterisation for the potential energy density, the relation ((18(1 becomes a 
differential equation for the Hubble parameter which can be integrated analytically. Thus determinations of H(z) 
can be used to constrain V{z) directly. Since effectively one will always be dealing with observations covering a 
finite redshift range, by an appropiate linear transformation in the redshift variable, we can always work in the 
interval [—1, 1], where we know the set of Chebyshev polynomials provide a complete orthonormal set of functions, 
i.e. any function in the interval can be expressed as a linear combination of Chebyshev polynomials. Moreover, these 
approximating polynomials have the smallest maximum deviation from the true function at any given order, and 
provide a well-defined estimate of the error introduced in the truncation of the expansion at a finite order. We point 
out, in passing, that such a parameterisation could be used in any other attempts considered in the literature where 
it was the equation of state the observable being parameterised by its redshift dependence. 

Using observations of passively evolving galaxies we obtain measurements of the Hubble parameter at 9 different 
redshifts. We use these determinations to constrain the first three coefficients in the Chebyshev expansion of the 
potential. For comparison we repeat the analysis using recent supernovae data, which give us an integral of the Hubble 
parameter. We find that the reconstructed potentials from both data sets are consistent, giving some confidence that 
the results are not heavily plagued by systematic errors. The standard LCDM model is consistent with current data 
at the lsigma level. We show that future data obtained from the Atacama Cosmology Telescope will greatly improve 
the constraints. 

Since a cosmological constant is a good fit to the observations we asked the question of how to generically describe 
small deviations from this simple scenario. It is clear that even if the nature of dark energy might be related to a 
dynamical field, distinguishing such a scenario from a real cosmological constant will be an extraordinary experimental 
challenge, as the dark energy potential can be arbitrarily close to a constant. 

We thus analysed the conditions for the dark energy field to "slow-roll" in the presence of matter, thus enabling a 
dynamical dark energy to get arbitrarily close to a cosmological constant ((4711 . By expanding the potential in Taylor 
series for A q (t) = q(t) — q(0) we derived the generalisation of the standard slow-roll conditions used in inflation, in 
the presence of matter, which translate into conditions that the functional V[q] must satisfy to explain the maximum 
deviation allowed from a cosmological constant. 

Even though in this paper we focused on single canonically normalised scalar field, it should be clear that it is 
straightforward to apply our formalism to an arbitrary number of them, not necessarily being canonically normalised. 
However, the observables quantities are the matter density, the Hubble parameter and its derivatives, which depend 
on the full kinetic and potential energy densities of the scalar field sector, and are insensitive to whether these values 
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are given by the superposition of more than one field. This opens up the question of whether one would be able to 
determine, experimentally, the existence of more than one rolling scalar field. In other words, to which extent it is 
possible to disentagle the full kinetic/potential energy of a superposition of scalar fields into the kinetic/potential 
energies of their components? 

We have illustrated the enormous experimental challenges of reconstructing V[q] from V(z) just for a single scalar 
field. The task is even harder for more than one field as there are more derivative directions to consider, dV/dqi, and 
direct experimental observables depend only on the time derivatives of the full potential energy, i.e.dV/dt = dV/dqi qi. 
Such a disentanglement seems extremely challenging, if not impossible, at least from the perspective of the formalism 
developed here. In this context we can say that our formalism enables one to reconstruct the properties (potential 
and kinetic energy), of an "effective" field. 

Appendix A 

Chebyshev polynomials can be computed using the recursion relation: 

T n+ i(x) = 2xT n (x) - T n -\(x) for n > 1 (73) 

where 

T (s) = l, T 1 (x) = x (74) 

Thus T n (x) will have the following structure: T n (x) = a x n + a 2 x n ^ 2 + a^x n ~ A + so for example: 

T 2 {x) = 2x 2 -I (75) 
T 3 (x) = 4a; 3 -3a; 
T 4 (x) = 8a; 4 - 8x 2 + 1 

(76) 

The integral on the RHS of Eq. (|28(l will involve a series of integrals of the type: 

/2z/z max -l 
x n (a + bx)- 7 dx (77) 

In particular, using the recursion relation of the Chebyshev polynomials we find that: 

T (x) = 1 F (x) = I 

Ti(x) =x Fx(x) = h ■ 

T 2 (x) = 2x 2 -1 => F 2 (x) = [2J 2 -7 ] • 

T 3 (x) = Ax 3 -3x => F 3 (x) = [47 3 - 3h] . 

Ti(x) = 8x 4 - 8x 2 + 1 => Fi(x) = [8/4-8/2 +io] • 

T 5 {x) = 16a; 5 - 20a; 3 + 5a; ^ F 5 (x) = [I6I5 - 20J 3 + 5/i] . 

T 6 (x) = 32a; 6 - 48x 4 + 18a; 2 - 1 => F 6 (x) = [327 6 - 48J 4 + 187 2 - J ] . 
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where 



, 2z 
L = 1 



l(l + z)- 6 (n-5)(l + z)- 5 (n-5)(n-4)(l + z)- 4 

(78j 



6 a 30a 2 120a 3 

(n - 5)(n - 4)(n - 3)(1 + z)" 3 (n - 5) • • • (n - 2)(1 + z)~ 2 (n - 5) • • • (n - 1)(1 + z) 



(-1) 



n+l 



360a 4 720a 5 720a 6 

1 (n-5) (n-5)(n-4) (ra - 5)(n - 4)(n. - 3) (n - 5) • • • (n - 2) 
6a ~ 30a 2 + 120a 3 360a 4 + 720a 5 



(n- 5) • • • (n- 1) 



720a 6 



720a 6 



and 



2z/z max -l x n 



1 



(a + bx) 



where 



40) = ~o n T \a ~ 1 

In the case of Eq. (|43() G n arc defined as: 

T {x) = 1 =► G (x) = Jo 



Tx{x) = x => Gi(x) = Ji . 
T 2 (x) = 2x 2 -1 => G 2 (x) = [2J 2 -J ]. 
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(79) 
(80) 



(n - 5) --- n J„ (81) 



(82) 



J„ = / 7 — dx (83) 



These integrals are given by 

■m—Q ^ ' 

but can also be obtained using the recursion relation 

J n = \ [(2z/z max - 1)" - (-1)"] - \ J„_! (85) 
no b 



J = \ log(l + z) (86) 
o 

Integrals /„ can also be obtained by recursive relation [^3 

j _ (l + z)- 6 (2z/z max - !)"-(-!)" n(l + w/2) 
(n - 6) j2 (n - 6)z max /2 

where 
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